Sharp signature of DDW quantum critical point in the Hall coefficient of the cuprates 
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We study the behavior of the Hall coeSicient, Rh, in a system exhibiting d^2_y2 density-wave 
(DDW) order in a regime in which the carrier concentration, x, is tuned to approach a quantum 
critical point at which the order is destroyed. At the mean-field level, we find that nnaii = I/Rh 
evinces a sharp signature of the transition. There is a kink in riHaii at the critical value of the 
carrier concentration, Xc; as the critical point is approached from the ordered side, the slope of nnaii 
diverges. Hall transport experiments in the cuprates, at high magnetic fields sufficient to destroy 
superconductivity, should reveal this effect. 



PACS numbers: 

There has recently been great interest in the possibil- 
ity that the "pseudogap" state of the cuprates is, in fact, 
an ordered broken-symmetry state . The interplay be- 
tween this ordered state and superconductivity is conjec- 
tured to lead to the suppression and eventual demise of 
superconductivity at low dopings. At a critical doping, 
Xc, which is near - but not necessarily precisely at - op- 
timal doping, the order associated with the pseudogap is 
expected to evaporate. In underdoped cuprates, the two 
orders are expected to coexist. 

An interesting way to establish the reality of this new 
order, which is hypothesized to be a (1^2 _y2 density- wave 
1^ (DDW), is to destroy the superconductivity in the 
nearby region of the phase diagram If the pseu- 
dogap is due to this order, it will survive, but if, on 
the other hand, the pseudogap were due to supercon- 
ducting fluctuations, it would collapse. The experiments 
by Boebinger and his collaborators in which supercon- 
ductivity is destroyed by 60 T pulsed fields are ideally 
suited for this purpose Q| . The effect of such magnetic 
fields, which are high enough to destroy superconductiv- 
ity, has been shown to have little effect on DDW order 
1^, ensuring conditions under which DDW order can be 
studied in its pristine form. In particular, the transi- 
tion between the DDW at lower carrier concentration to 
the "normal" state at higher carrier concentration can be 
studied at asymptotically low temperatures. An impor- 
tant question is whether or not this transition leaves a 
sharp signature in a measurable property, with emphasis 
on the word sharp. Boebinger and his collaborators ||] 
have addressed precisely this question by measuring the 
low-temperature Hall constant Rh or, rather, its inverse, 
n-iiedi = ^/Rh = (o'xx)'^ / o'xy In conventional Drude the- 
ory, it is equal to ne/B, where n is the carrier density, e is 
the electron charge, and B is the applied magnetic field. 
Preliminary measurements indicate indeed a sharp 
signature in rtHaii at a carrier concentration Xc 0.15. 
While these difficult and important experiments need to 
be firmly established, we wish to discuss this question 
theoretically from the DDW perspective. Thus, exper- 
imental results provide the context and the motivation 



for our theoretical work. 

We study the Hall number, nnaii, close to the quan- 
tum phase transition at zero temperature. The theo- 
retical framework is an effective Hartree-Fock Hamilto- 
nian that captures the broken symmetry of the DDW 
order. The contribution of nodal quasiparticles, arising 
from this Hamiltonian, to electrical and thermal trans- 
port was previously studied, and it was shown how they 
can be used to detect DDW order . However, near the 
quantum critical point at x — Xc, where DDW order dis- 
appears, the DDW gap is small, and it is not sufficient 
to focus solely on the nodal quasiparticles; quasiparti- 
cles far from the nodes are equally important. In this 
paper, we report the results of a calculation which takes 
into account all of the relevant parts of the Fermi surface. 
We consider the Hall coefficient in linear response theory, 
using the Kubo formulae. We find that the Kubo formu- 
lae lead to the same results as the Boltzmann equation 
in the weak scattering limit, unlike in the case of nodal 
quasiparticles, where the Boltzmann equation cannot be 
naively used. Consistent with the Hartree-Fock effective 
Hamiltonian, we assume a mean-field dependence of the 
magnitude of the zero-temperature DDW gap, A(a;), on 
Xc — x, A{x) = Ao{xc — xY^^ , for x < Xc', here on, we will 
drop the argument of A(a;), if there is no danger of con- 
fusion. We find a kink in riHaii against x. As x —> x~ , the 
slope of riuaiiix) diverges; for a; — > a;+, the slope remains 
finite. As we discuss below, the change in the slope is due 
to the contribution of the "hot spots" , where the Fermi 
surface crosses the reduced Brillouin zone (or magnetic 
zone) boundary ||^. 

It is clear that fluctuation effects can not be neglected 
at the quantum critical point On dimensional 

grounds, however, there is some reason for believing that 
the quantum critical region is very narrow and is of or- 
der {Aq / Ep)'^ , where Ep is the Fermi energy. It would 
indeed be very difficult to explore the critical region, as 
it would require preparing samples with an extraordi- 
nary degree of control of doping. Thus, the choice of 
the mean-field critical exponent is eminently reasonable. 
As an ad hoc procedure, we have also discussed other 
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possible dependences of A on Xc — x, only because such 
choices have appeared in the literature jo). We want to 
warn the reader, however, that these choices are merely 
empirical and no theoretical justifications can be offered. 
The actual theoretical problem of incorporating fluctu- 
ations at the quantum critical point is more complex, 
which is squarely beyond the scope of the present paper; 
we further remark on this issue later. 

First, we use the Kubo formula to derive the conduc- 
tivity and Hall conductivity of a system of electrons with 
DDW order, assuming that the only source of scattering 
is impurities. 'Residual' interactions which remain after 
the development of DDW order are neglected. We begin 
with the mean field Hamiltonian for the DDW state: 

H ^[{f^i, - M) cL^ka + «Ak4„Ck+Qa + h. c] (1) 
k,a 

where Ck is the annihilation operator for an electron of 
spin a in the z-direction and momentum k. The sin- 
gle particle spectrum on the square lattice with nearest- 
neighbor hopping t and next-neighbor hopping t' is ek = 
4t' cos kx cos ky — 2t{coskx + cosky). Ak = A(cosfc2: — 
cos ky) is the d- wave order parameter of DDW state 
and the vector Q = (tTjTt). We have set the lat- 
tice spacing to be unity. We can express the Hamilto- 
nian in terms of a 2-component quasiparticle operator: 
a ~ ('-ka' ~*''k-i-Qa)' ^'^'^ then diagonalize this simple 
2x2 Hamiltonian, to get 



k,CK 











Xk,Q, 



(2) 



here. 



eik 



where E±{k) = e2k ± 
—2t{coskx + cosky), e2k = At' cos kx cos ky . The two- 
component quasiparticle operator Xka is a unitarily re- 
lated to 5'kQ, and the sum is over the reduced Brilloun 
zone (RBZ). 

Since the field x is the superposition of two charge — e 
fields, it, too, is a charge — e field. Hence, the electrical 
current operators are given by 



'HxL 

k.a 



dE+{k) 

, y 






Xka 



(3) 



In the Kubo approach, conductivities are obtained by 
applying weak electrical and magnetic fields to the sys- 
tem: E = EQXCOs{ujt), B — qzAQsm{qy). At the linear 
response level, the conductivity is given by: 



o-xx = — Imn2(iw„ 
ui 



iS,q = Q,T) 



(4) 



where Hz (ic^„,q,r) = /j'dre'-"- < T,jx(q, r)j,(q, r) > 
is the Fourier transform of the imaginary time-ordered 
current-current correlation function, and aj„ is the Mat- 
subara frequency (/3 is the inverse temperature). The 
Hall conductivity is given by pXl]: 



axy{uj, T) = lim — ReH3(icj„ uj + iS, qy, T), (5) 
9— »o Loq 

where 

H3(zc^„,q,T) = 

drdr'e'^-^ < Trjyiq, r)j,(0, 0)j,(-q, r') > (6) 

/o 

In the limit of large scattering time, Ts, we make the 
usual approximations [ pH , and taking g — + 0, w — > 
limits, we get 



(Tx,i = e^Bri 



d^k 
(2^ 

dE+{k )\^ d^E+{k) 
dkx 



dE+[k) dE+{k) d^E+{k) 



dkr 



dkl 



d{E+{k)-fi) + iE+^E^), (7) 



and 



e r, 



d^k fdE+{ky^ 



(27r)2 

{E+^E^] 



dkx 



8{E+{k)-p) 



(8) 



These are precisely the same formulae that are obtained 
from the Boltzmann equation |p^ . 

With formulae (|^) and (||) in hand, we can proceed to 
calculate the Hall coefficient. Above the critical hole con- 
centration [xc ^ 0.2), there is no DDW order, and riHaii 
is obtained from the same formulae, but with the quasi- 
particle dispersions E± (k) replaced by the tight-binding 
electron band structure ek. In this case, the sign of the 
Hall number is determined by the curvature of the tight 
binding electron Fermi surface. If the Fermi surface is 
closed around the center of the first Brillouin zone (the 
F-point), the Hall effect will be electron-like, and nnaii 
will be negative. On the other hand, if the Fermi surface 
closes around the corners of the zone (the Af-points), 
?^Haii is positive, and the response is hole-like. Band 
structure calculations ]l3[ |, corroborated by a wide spec- 
trum of phenomenology, indicate that the latter scenario 
occurs in the cuprates. This is a constraining factor in 
our choice of parmeters. 

Insofar as the behavior close to the critical point is 
concerned, the variation of is clearly smoother than the 
variation of A, and we can treat ^ as constant, thereby 
neglecting the smooth variation of the Hall coefficient 
which is generated by the variation of /i with x. The re- 
sults shown below correspond to t = 0.3 eV, t' /t — 0.3, 
/i = — i, and Xc = 0.2. The precise value of Xc is unim- 
portant in our calculation. Its value is non-universal and 
could vary between materials. For example, it is about 
0.15 in Ref. §, while it is about 0.2 in Ref. §]. Because 
we hold constant around Xc, we will obtain a Hall co- 
efficient, which is also a constant for x > Xc- Later, we 
will comment on the correction to this assumption. The 
order parameter. A, is assumed to scale with a mean field 
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exponent near the critical point: A(a;) = Ao{xc — xf-^"^ ■ 
For a representative value, we have chosen Aq = 0.03 
eV. We wish to emphasize that the results of our calcula- 
tion are robust with respect to the choice of parameters, 
as long as the considerations pertaining to the hot spots 
described below hold. We have explicitly tested this by 
choosing a wide range of parameters. 

Let us now analyze equations and (||) close to the 
critical point. The formulae involve integrals over the 
Fermi surface of first and second derivatives of the quasi- 
particle dispersion. Let us first consider a typical deriva- 
tive, such as ; which appears in the integrand on the 
ordered side of Xc, for small values of — xd. 



dE k 

— - — = — 4t' sin fcr cos fc„ 



. f (cos 
sm kx 

[4i^(cos kx -I- cos kyY 



kx — cos ky) -f- 4i^(cos kx -I- cos ky) 

^ ' - Ai(cos fca; — cos fcj,)2]i 



(9) 



For small values of \x — Xf\ on the ordered side, A is 
negligible over most of the region of integration. How- 
ever, this is not true at those points on the Fermi sur- 
face where t(coskx + cos ky) is even smaller than A; 
at these points, A is important. This will be the case 
at the points where the Fermi surface crosses the RBZ 
boundary - the "hot spots" |^ - and finite regions 
\coskx + cos fcj^ I / 1 cos fca; — COS fcj^ | < A /At about them. 
These regions of integration will give rise to a slope dis- 
continuity in riHaii- 

To illustrate this, we consider the quantity along 
the boundary of the RBZ: 



dkx 



= —At' sin kx cos ky 



A 



■ sin kr 



(10) 



On the other hand, in the disordered phase for x > Xc 



A is identically zero, so 



9^4 



takes the following value 



along the boundary of the reduced Brillouin zone: 
dE 

— = —At' sin kx cos k,, — 2t sin kr 

dkx X y X 



(11) 



Thus, in the A limit, this quantity has a discontinu- 
ity of order t at the RBZ edges. 
Consider, now, ^^^g^ for A ^ 0: 



dA 



2u^ 



dA 



— . We would like this quantity to zeroeth order in 

A. Thus, we can take A = in the factors and 
Inside the derivatives, we only need to keep the linear in 
A terms in axx and axy Hence, the computation reduces 
to understanding the leading A dependence of axx and 
axy Consider the integral of Eq. the integral of Eq. ^ 
is similar. We can calculate (^) for A small as the A = 
value plus a correction term coming from the integral 
over the region around the hot spot. From ( pO| ) and (11), 
we see that the integrand - which is a product of three 
such derivatives - is finite in the A = limit. Hence, the 



leading A dependence comes simply from the size of the 
integration region, \coskx + cos fcj, | / 1 cos fc^, — cos ky\ < 
A /At, which is tA. 

Thus, we find that the leading A dependence as x — > 
x~ is linear. As A ^ 0, this term vanishes, so the Hall 
constant is continuous. However, its derivative with re- 
spect to X will not vanish The slope of nnaii 
(normalized to its value at Xc) will actually diverge for 
X ^ x^ . In Fig. 1^, we have plotted the calculated val- 
ues of TiHaii with X close to Xc on either sides, with the 
mean- field dependence of A (a;). As we see, the extra 
contribution resulting from non-zero A is negative. 
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FIG. 1. The Hall number calculated as a function of dop- 
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at the quantum critical 



ing X. Note that 
point as a; ^ x~ . nnaii was normalized to its value at x^- 

The effect is further illuminated by considering the evo- 
lution of the Fermi surface itself. As the system crosses 
from the disordered side of the transition to the DDW 
side, the Fermi surface becomes disconnected precisely 
at the hot spots. In the extreme case, t' = fi — 0, where 
the tight-binding Fermi surface is the RBZ boundary, 
this occurs everywhere on the Fermi surface. The Fermi 
surface is reduced to the four nodal Fermi points even 
for infinitesimal A (it is disconnected everywhere) . How- 
ever, for any finite fi,t', the disconnection takes place 
only at the eight crossing points, the hot spots. The 
typical DDW band structure is shown in Fig. ||. 

The nature of the discontinuity in the slope of /iHaii de- 
pends on the critical behavior of the DDW single-particle 
gap. One may wonder, what would happen if the gap be- 
haved more generally with a critical exponent /?' given by 

A{x)^Ao{x,-xf7 (12) 

Superficially, the slope will vary as \x — Xc\^ as the 
transition is approached from the ordered side, if we plug 
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FIG. 2. The DDW bandstructure plotted along special 
directions of the Brill ouin zone of the square lattice; the 
point group symmetry is 4mm. The parmeters a re described 
in the text. The upper band is the dashed curve and the 
lower band is the dashed-dotted curve; the chemical potential 
is noted by a solid line. Note that the chemical potential lies 
in the upper band at the X point. 



this A(x) in our equations. If /?' < 1, as is the case in 
mean- field theory (/3' = 1/2), then the slope will diverge 
as the transition is approached from the ordered side. 
If /?' = 1, as Loram et al. suggests then there is a 
finite slope discontinuity. If /3' > 1, then the slope is con- 
tinuous at the transition and non-analyticity only show^s 
up in higher derivatives of riHaii with respect to x. It 
is inconsistent, however, to use a non-mean field expo- 
nent in the context of the Hartree-Fock Hamiltonian, as 
the exponent f3' ^ 1/2 must arise from fiuctuations at 
the quantum critical point for which we can no longer 
meaningfully use the Hartree-Fock Hamiltonian, or its 
quasiparticle excitations. As remarked earlier, the criti- 
cal region is expected to be narrow and not of relevance, 
given the precision of present experiments. 

To go consistently beyond mean-field theory, we must 
not only take the correct value of /3 - which may be 
the 3D Ising exponent since the DDW order parameter 
breaks a Z2 symmetry - but also consider the scattering 
of quasiparticles by critical fiuctuations. 

Although our calculation was performed for DDW or- 
der, superficially, identical behavior would be obtained 
in the case of any order parameter at Q = (7r,7r): there 
would be linear in A contributions originating at the hot 
spots. The important task in this case is to phenomeno- 
logically motivate such a order parameter transition. We 
find no compelling evidence for the existence of other 
order parameters, such as a commensurate spin-density 
wave, or a triplet d-density wave setting in at a; = Xc. Hy- 
pothetical order parameters at other wavevectors could 



have some signature, but it might not be linear in A, 
so these order parameters may require smaller values of 
/3' in order to manifest themselves as kinks in riHaii or 
else they may only show up in higher derivatives of riuaii ■ 
So, it is crucial to precisely determine the experimental 
nature of the signature at Xc in riHaii- 

In a real experiment, the Hall coefficient will be mea- 
sured at a finite temperature. Its value will depend on 
the order of the limits A ^ 0, T ^ 0. As long as A > T, 
the result of the previous section will be observed. How- 
ever, if it is possible to approach very close to Xc, so 
that A < T, then A is replaced by A'^/T. If /?' = 1/2, 
a slope discontinuity remains. For /?' > 1/2, the tran- 
sition is rounded. To summarize, an experiment at low 
non-zero temperature will observe a rounding of the tran- 
sition within a very small window Xc — x < T^lf^ . 

In our calculation, we have ignored the variation of fi 
with X. If we allow /i to vary with x in the manner dic- 
tated by the single-particle spectrum, then the result of 
Fig. |l| is superposed on a large but finite sloping back- 
ground. In this case, rtHaii will not be a constant for 
X > Xc- However, there are some indications that the 
chemical potential varies more slowly in the underdoped 
regime than one might expect [Q. At any rate, this 
background is non-universal and determined by separate 
physics from that which governs the critical properties. 
It is really a separate issue, and we do not attempt to 
model it here. 
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